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Abstract: We examine a possibility to explain the signals of dark matter direct detection 
from DAMA and CoGeNT experiments, under the framework of supersymmetric extension 
of the standard model. For this purpose, we introduce four Higgs doublet fields and assume 
that one pair of Higgs fields which have very small vacuum expectation values gives sizable 
effects for annihilation and scattering processes of dark matter. We show that the preferred 
parameter regions for DAMA and CoGeNT results can be simultaneously explained by 
supersymmetric four-Higgs doublet model with the parameters consistent with the observed 
value of dark matter relic abundance. The extra Higgs fields introduced as an explanation 
of the light dark matter scenario also explain the Wjj anomaly reported by CDF. 
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1 Introduction 

The standard model is very successful in describing the interactions of elementary parti- 
cles. However, from a cosmological point of view, we have not reached the theory which 
can explain the whole nature of the universe. In particular, a number of observations 
suggest that most of the mass in the universe is composed of non-baryonic dark matter. 
Interestingly, some experiments claim that they observed a signal compatible with the 
light dark matter candidate by the direct detection search. Among them, DAMA [1] and 
CoGeNT [2] experiments have reported an annual modulation signals which support the 
existence of dark matter. Moreover, CRESST experiment recently reported the observation 
of the signal events with more than 4(T significance [3]. Although other experiments such 
as XENONIOO [4] and CDMS-II [5] have already reported the null result which exclude 
the preferred region of the DAMA and CoGeNT results, we put emphasis on the fact that 
there is 8.9a signature of an annual modulation by DAMA experiment [6]. For the validity 
of the analysis of DAMA experiment, some discussions are found [7] [8] [9] and it should be 
investigated more carefully in the future analysis. 

Our goal in this paper is to construct a model which is capable of explaining the light 
dark matter signals claimed by DAMA and CoGeNT experiments. There are a number of 
works which investigate the model to explain these results [10] [11] [12] [13] [14] [15] [16] [17] [18], 
and we try to construct a concrete model as an extension of the standard model. For this 
purpose, we need to introduce a stable neutral field which has a coupling large enough 
to generate the observed rate of dark matter direct detection. In other words, the key 
ingredient for the model building is a coupling between the dark matter candidate particle 
and the field which mediate the dark matter-nucleus scattering process. 
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As a model to realize a dark matter candidate, the minimal supersymmetric standard 
model (MSSM) is well motivated because the lightest supersymmetric particle (LSP) is sta- 
ble by virtue of i?-parity. In fact, there are some attempts to interpret the DAM A result 
under the framework of the MSSM [19] [20] [21]. According to these study, large Yukawa 
coupling and light pseudoscalar which mediate for the process of neutralino pair annihila- 
tion are required to obtain the realistic value of neutralino relic abundance. However, such 
scenario with large tan /3 and small pseudoscalar mass Ma is severely constrained by the 
direct searches at the Tevatron and the LHC [22] [23] [24] [25]. 

With this knowledge, we consider a model of the extension of the MSSM which has 
an extra pair of Higgs doublet superfields. There are already some works which introduce 
multi-Higgs fields as an extension of the MSSM [26] [27] [28] [29] [30] [31] [32] and we try to 
construct a model of light dark matter under the framework of supersymmetric four-Higgs 
doublet model. That is, if the LSP neutralino include some contribution from the extra 
higgsino components, this lightest neutralino couples to the scalar components of the extra 
Higgs fields. Therefore, we can obtain a sizable contribution for the annihilation and 
scattering process of dark matter and quarks if there are large Yukawa couplings for the 
quarks and extra Higgs fields. Note that we can achieve large Yukawa couplings for extra 
Higgs fields assuming that these fields obtain quite small vacuum expectation values (VEVs) 
and evade the usual Higgs search constraints. 

Although it is not difficult to construct a model which is suitable for the favored param- 
eter region of DAM A experiment, reconciling the contradiction with the null experiments 
such as XENONIOO and CDMS-II is problematic. Moreover, there is some discrepancy 
between the favored parameter regions of DAMA and CoGeNT experiments. Since the 
favored region for DAMA result is sensitive to the quenching factor of sodium QNa; there 
are some arguments for the choice of this quenching factor. If we allow larger value of 
quenching factor than the default value Qnu = 0.3 it 0.03, DAMA favored region extend 
to include the smaller mass range which relax the discrepancy with the CoGeNT result 
[33] [34]. 

As the possible solution for the contradiction among the experiments, there exist some 
works which propose an isospin violation in the dark matter-nucleus cross section [17] [18]. 
In fact, taking the ratio of neutron to proton couplings A„/Ap ~ —0.7 and requiring large 
sodium quenching factor QNaj DAMA and CoGeNT preferred region and the null result 
from XENONIOO become consistent for a dark matter with mass ~ 8 GeV [33] [34]. 
In our model, we can easily achieve this isospin violation by adjusting the extra Yukawa 
couplings for up-quark and down-quark as suggested in [18]. Since it is still impossible 
to reconcile these results with the limit from the CDMS experiment, we cannot claim the 
success for the explanation of all experiments. Even so, it is interesting that we can explain 
other experiments by the simple extension of the MSSM, and the future experiments should 
test the model by settling the contradiction. 

In our model, we should have extra Higgs fields with mass of 0(100 GeV) to acquire 
the dark matter-nucleus cross section required to explain the DAMA and CoGeNT results. 
Therefore, it may be possible to detect this extra Higgs field produced by the collider 
experiments. Interestingly, CDF Collaboration have reported an data on the dijet mass 
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distribution of the W + jj channel which indicate a 3.2a" excess around 150 GeV [35]. 
One way to interpret this excess is to introduce new particles which couple to quarks 
at tree-level [36] [37] and the extra Higgs fields we introduce might be a solution for this 
anomaly. Although this anomaly is somewhat questionable because it is not confirmed 
by DO experiment, the connection with the dark matter search predicted by our model is 
interesting and should be tested by LHC. 

The outline of this paper is as follows. In Section 2 we introduce the basic setup 
of a model of dark matter which is based on supersymmetric four-Higgs doublet model. 
In Section 3 we consider the thermal relic abundance of dark matter in our model with 
the neutralino LSP scenario and compute the direct detection rate of dark matter. In 
Section 4 we investigate the constraints and discovery potential from other experiments. 
Our conclusion is given in Section 5. 

2 Supersymmetric four-Higgs doublet model 

Supersymmetric four-Higgs doublet model is a simple extension of the MSSM by extra 
Higgs doublet superfields which have the same quantum numbers as the original up-type 
Higgs Hu and down- type Higgs Hd- In this section, we illustrate the detail of the model 
which satisfy the requirements for explaining the properties of dark matter. 

To describe Higgs fields in this model, we use a following notation for two pairs of 
up-type and down-type Higgs superfields: 



where we assume that the masses of quarks and leptons come mainly from the VEVs of 
Hui and Hdi- In other words, we introduce Hu2 and Hd2 as fields with very small VEVs 
compared with Hui and H^i- From now on, we will refer to the Higgs doublets Hui, H^i 
as "original Higgs" which correspond to the fields originally contained in the MSSM and 
the additional Higgs doublets Hu2 , H(i2 as "extra Higgs" . 

Now we can write the superpotential terms involving these fields as 



where a, b (= 1,2,3) are fiavor indices for the matter superfields. Other terms such 
as L°^Hui are assumed to be forbidden by i?-parity. Once Higgs fields obtain the VEVs 
which break the electroweak symmetry, each field contributes to the masses of quarks and 
leptons. Since we assume that Hu2 and Hd2 obtain very small VEVs (of order less than a 
few MeV), they do not give the large contribution to the fermion masses even if we assign 
0(1) couplings for the Yukawa interactions among the matter fields and extra Higgs fields. 
This assumption is crucial to explain the results of dark matter search as we will show in 
the next section. 

The soft supersymmetry breaking terms which are relevant to the Higgs scalar potential 

are 




(2.1) 




1=1,2 




(2.3) 



i=l,2 



1=1,2 
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Then the scalar potential of Higgs fields can be written as 



j=l,2 



+ [(^11^21 + fJ-u 1^-22) {H^IH ^2 + H^*H^2) 

+ (/^11/"12 + ^^lA^22)«*i/,°2 + ^d7^d2) + C.C] 



+ 



+ 



i=l,2 



^2 r 



+ 



di) 



|2 

I 



(2.4) 



Although the couplings /Xjj and bij are complex parameters in general, from now on we 
take these parameters real for the sake of simplicity. 

As was pointed out in [38] , the VEVs of Higgs fields in this model is always real when 
all of the parameters in the Higgs potential are real at tree-level. Now we assume that 
neutral Higgs fields obtain real VEVs = Vui / 0, {H^^ = Vdi / and the VEVs of 

charged Higgs fields vanish so as not to break the electromagnetic symmetry. So we can 
parameterize these VEVs as follows: 

Vui = visin(3i, Vdi = vi cos f3i, f ^2 = ^1 tan a; sin /32 , ^^2 = ^^i tan a; cos /32 (2-5) 

where v'^ = vf{l + tan^w) = 2M^/(gf^ + g''^). Since we consider the situation where the 
extra Higgs fields obtain very small VEVs, we put tan a; <C 1. 

Demanding that the first derivatives of the Higgs potential with respect to the neutral 
Higgs vanish, dV/dH^^ = dV/dH^^ = 0, we obtain 

= /^ii + /^12 +"^ni + (/^ii/"2i + ^12/^22) tana; sin~^ /3i sin/32 
— 611 tan~"^ /?! — 612 tantj sin^"*^ /3i cos (32 



+ 



M2 



[sin^ /?! — cos^ /3i + tan^ w(sin^ (32 — cos^ /32) 



2(1 + tan2 w) 

= /^ii + /^2i + (/"11M12 + /U2i/i22)tana;cos"^/3i cos/32 

— 611 tan /3i — 621 tanw cos~^ Pi sin (32 



(2.6) 



[O O 0/0 On 

sin /?! — cos /3i + tan 6t;(sin /32 — cos (32) 

= /^2i + fJ'22 + "^n2 + (/^ii/«2i + ^12/^22) tan~^ w sin /3i sin"-*^ (^2 
— 622 tan~"^ /32 — 621 tan""*^ w cos (3i sin~^ /32 



(2.7) 



/Vf2 

7 rO O 0/0 On 

H — ; ;5 — r sill /3i — cos /3i + tan 6t;(sin /32 — cos /32) 

2(1 + tan"= CO) '- 

= /^12 + A*22 + ^d2 + (/^ii/"i2 + /U21M22) tan"^ w cos /3i cos"^ (32 



(2i 



-4- 



(2.9) 



— 622 tan 1^2 — bi2 tan ^ uj sin /3i cos ^ /32 

- -7- ^^^s — r [sin^ /?! - cos^ /3i + tan^ w(sin^ /32 - cos^ /32)1 . 

2(1 + tan^ oj) 

For tanu <^ 1, equations (2.6) and (2.7) reduce to the approximate constraints 

~ /xfi + /x?2 + mil - bii tan-^ /^i " ^ ^os 2/3i, (2.10) 



/Vf2 

— /^ii + /"21 + "'-di ~ ^11 tan/3i + cos2/3i, 



(2.11) 



and to satisfy the constraints (2.8) and (2.9), there are four possible choices of the param- 



eters, I.e., 



1. 612, 621 ~ and //II, /U22 — 0. 

2. 612, 621 — and /U12, At2i — 0. 

3. 612, 621 ^ and fiu ^ ±/i22, W2 ^ T/"2i- 

4. 612 ^ (/iiim2 + //2i^22) tan~^ /3i and 621 ^ (/^ii//2i + /t/12/^22) tan/3i, 

because of large tan^^ tj. Since we need a tuning of the parameters for the case 3 and 4, 
we will adopt the case 1 and 2 for the subsequent analysis. 

Now let us consider the masses of the particles related with the Higgs sector. Since 
neutralinos contain the higgsino components corresponding to the extra Higgs, we write 
the masses of neutralinos in the basis = {B, W^, H^i, H^i,H^2^ ^u2) 

( Ml -cpswMz spswMz \ 

M2 cpcwMz -s^cwMz 

-cpswMz cpcwMz -^11 -^21 

spswMz -spcy/Mz -/Uii 
-1112 

V -^21 



-fJ-12 

-fJ-22 
-/^22 / 



(2.12) 



where we adopt the notation that 

s^ = sin/3i, c^ = cos/3i, 



ts = tan/3i 



(2.13) 



and neglect the elements which are proportional to tanw. Similarly, the mass matrix of 
charginos in the basis 4^ = {W~^ , H^^, H^2^ -^di' ^ds) 



M. 



c 



'0 

M 



M ~ 



/ M2 ^^2s,iMw 
\/2cpMw /ill ^J'2l 

\ /il2 ^X22, 



(2.14) 



and the spectrum of these particles can be modified from that of the MSSM if /ii2 and /i2i 
are not so small. 
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Next let us consider the Higgs scalar fields. For the real parts of the neutral Higgs 
fields {ReH^i,ReH^^,ReH^2^^^^d2) ^ have CP-even Higgs mass matrix 



Mi 



( but/ + Mis} 
-bii - Mlspcp 


V 



-611- 

biitp 











~b22 



\ 


-b22 

Ad J 



(2.15) 



and for the imaginary parts of the neutral fields (ImfZ"!]]^, Im//^^, Imff^2) I™-f^d2)' CP-odd 
Higgs mass matrix is 



(biitp 


611 





0\ 


bu 


but/ 














A„ 


^22 


V 





&22 


Ad/ 



(2.16) 



where A„ and A^ are defined as 
A„ 



f4l + M22 + "i«2 



C2/3, Ad = + ^22 + "^d2 + ^ 



C2/3- 



(2.17) 



Here we used (2.10) and (2.11) to simplify the matrix elements, and omitted the elements 
which are proportional to 612, 621 and tana;. It should be pointed out that these mass 
matrices are block diagonal and the structure of the blocks corresponding to the original 
Higgs H^^ and H^^ are similar to that of the MSSM. Moreover, for the blocks corresponding 
to the extra Higgs fields, the mass eigenvalues are degenerate between CP-even Higgs and 
CP-odd Higgs. So we find the mass eigenstates of the neutral Higgs fields; 



^u2 
^d2 




(2.18) 
(2.19) 



where Sa = sin a and = cos a represent the mixing of the real components of the original 
Higgs and = sin 7 and = cos 7 are the mixing angle of the extra Higgs. Here we define 
the mixing angle 7 so that the eigenstates corresponding to the lighter mass eigenvalue are 
</>g-^ and (pQi- Note that the sign of the mixing angles of the extra Higgs are opposite 
between CP-even components and CP-odd components. 

The mass matrix of charged Higgs fields [H^^ , H^^, ^ ^da) ^^"^ written in terms 
of that of CP-odd Higgs fields as 



Mi 



Ml + M,2 



w 




\ 



\ 



C2/3 

-C2pJ 



(2.20) 
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3 Neutralino dark matter scenario 



In this section, we investigate the very light (~ 10 GeV) neutrahno dark matter scenario 
under the framework of supersymmetric four-Higgs doublet model. We expect that this 
light neutralino explain the signals of dark matter by means of the exchange of extra Higgs 
fields for annihilation and scattering processes of neutralino. Especially, our model have an 
advantage that we can realize an isospin violation in the dark matter-nucleus cross section 
which is required to satisfy both of the signals of DAMA and CoGeNT [17] [18]. 

To discuss the neutralino dark matter scenario in this model, we first investigate the 
mixing structure of the lightest neutralino Xi- This lightest neutralino is guaranteed to be 
stable by virtue of ii-parity as long as this particle is the LSP. Let us write the lightest 
neutralino in terms of its components fields, bino B, wino and higgsinos i/^j, Hui- 

X? = A^ii^ + Ni2W^ + NnHdi + NuHui + A^i5i?d2 + A^i6^n2 (3.1) 

where is a unitary matrix which diagonalize the neutralino mass matrix, N*M^N~^ = 
Mp^. If we denote the mass of the lightest neutralino as m^, the coefficients Nu (i = 
1, . . . , 6) should satisfy the relation 

= (Ml - m^)Nu - c^swMzNi^ + spSwMzNi^, (3.2) 

= (Ms - m^)Ni2 + cpcwMzNi^i - spcwMzN^, (3.3) 

= -cpswMzNii + cpcwMzNi2 - m^Ni^ - /UnA^u - /UaiA^ie, (3.4) 

= spswMzNii - si3CwMzNi2 - HuNis - m^Nu - /Ui2A^i5, (3.5) 

= -/ii2A'i4 - m^A"i5 - ^22A^i6, (3.6) 

= -//2iA'i3 - ^22A^i5 - m^NiQ. (3.7) 

as in the case of the MSSM [39]. 

Meanwhile, we are interested in the situation where (1) /in, //22 — or (2) fii2, fJ-21 — 0. 
In the case 2, it is obvious that the lightest neutralino can contain very small components 
of the extra higgsino -^«2- Therefore, we adopt the case 1 and assume //n, fi22 

Ml <^ M^, /ii2, /X2i5 M2. The assumption Mi <C Mz is to realize the situation of the very 
light neutralino dark matter of 0(10 GeV). Since the masses of charginos are constrained 
by LEP II bound to be heavier than about 100 GeV, the lightest neutralino cannot be 
wino-like or higgsino-like. With this assumption, we can solve (3.2)-(3.7) to obtain 

Ni2 - ^^^^twNii « 1, (3.8) 
Ni3 ^ ^ATie < 1, (3.9) 

/^21 
TTl 

Nu ^ ^Ni5 « 1, (3.10) 



My 

Ni5^S0Sw—Nu, (3.11) 

My 

~ -cpsw—Nn, (3.12) 
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and 



lc^Nf,+Nf, + Nf, 




(3.13) 



The mass of the hghtest neutralino can be written as 




) 



-1 



■X — 



Ml 1 + 



(3.14) 



3.1 Dark matter relic abundance 

For a stable neutralino, we can estimate the neutralino thermal relic abundance by cal- 
culating the neutralino pair annihilation cross section. In this model, it is reasonable to 
assume that the extra Higgs exchange contribution is dominant for the process Xi Xi ~^ / / 
if the Yukawa couplings for the extra Higgs given in (2.2) are large enough. Especially, we 
will investigate the situation where the extra Higgs which correspond to the lighter mass 
eigenvalue {(pEi and (poi) dominate the contribution of the pair annihilation, and assume 
that the other eigenstates {(j)E2 and (j)o2) are heavy enough to be neglected. Since CP-even 
extra Higgs (pEi and CP-odd extra Higgs (poi are degenerate as pointed out in the previous 
section, we denote their masses as from now on. 

Although the extra Higgs exchange contribution for the neutralino pair annihilation 
process is generated from the Yukawa interactions in (2.2), such contributions for the flavor 
interactions are severely constrained by plenty of experiments. First, if there is a large 
Yukawa coupling for the electron and extra Higgs, it leads to the process e"*" e~ — )• cp^^^ and 
the constraint > 200 GeV comes from the LEP II experiment. Moreover, off-diagonal 
components of the Yukawa matrix for the quarks and extra Higgs are prohibited because 
such terms induce the flavor changing neutral current at tree-level. For the subsequent 
analysis, we simply assume that all Yukawa matrices for the quark sector are proportional 
to a unit matrix in a basis of quark mass eigenstates; 



VuL^ui {Hui) VuR ^ diag (m„, rric, mt), V^^Yai {Hai) VdR ~ diag (m^, ms,mh) (3.16) 

and we parameterize the size of extra Yukawa couplings by y'^ and y'^. To avoid the 
constraint from the extra Higgs production process by LEP II, we simply assume that 
there is no Yukawa interaction term between leptons and extra Higgs. 

With this assumption, we calculate the process Xi Xi ~^ / / where / represents quark 
fields which contribute to the annihilation process. Note that we neglect the quark masses 
and consider the contribution which comes from the exchange of extra Higgs fields cpEi and 
(poi- For the process with the final state up-type quarks, Xi Xi ^ uu, we can compute the 
thermally averaged annihilation cross section (crv)^ which is given by 



(3.15) 



where 




(3.17) 
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Here, we leave only an s-wave annihilation contribution because p-wave contribution is 
suppressed by a factor = T/m^. Similarly, for the process involving down-type quarks, 
XiXi ^ d d, we obtain 

i<^^)d ^ ^^9'yiNf,s^^{N,,s, - N^^c^f. (3.18) 

Therefore they sum up the result 

(av) = Nu ■ (av)^ + iV, • {av)^ = —^g'^Nf,{2y'^c'^ + 3y2s'^){Ni,s^ - Niec^f (3.19) 

where we take = 2 and = 3 because quark fields other than top quark give the con- 
tribution for this annihilation process. Since the annihilation cross section is proportional 
to m^/m^, this cross section tends to be small for the very light neutralino. 

Using the annihilation cross section, the relic abundance of the neutralino can be 
obtained as [40] 



2 ^, Xf ■ 1.07 X 109(GeV-^ 

"■72 '/ r 

gj mpi{av} 



^xh' ^ ' \„ ^ (3.20) 



where x/ is a freeze-out temperature and is an effective massless degree of freedom. Now 
we can determine y'^ from the experimental value Q.^h? ~ 0.1. For the case \y'^\ > 
required value of can be estimated to be 

where we use the values (7* ~ 8 and Xf 20. Note that this Yukawa coupling tends to 
blow up above the mass scale due to the contribution to the renormalization group 
equation from the loop mediated by the extra Higgs. So this imply that there might be 
some other new physics around a few TeV. But we will not discuss further about this topic 
in this paper. 

3.2 Direct detection cross section 

As we mentioned in the introduction, direct detection experiment is interesting because 
DAMA and CoGeNT experiments have reported that they observed the annual modulation 
of the signals which is compatible with a light dark matter scenario. So we investigate 
whether our model can explain these signals by computing the rate of the scattering between 
dark matter and nucleus. Requiring that the proper value of the relic abundance of dark 
matter ^-^Ji^ ~ 0.1 is generated, we can expect to derive some prediction for the direct 
detection rate. For the explanation of the results of DAMA and CoGeNT, an isospin 
violation plays the crucial role as we show later. 

To calculate the detection rate of neutralino-nucleus elastic scattering, let us consider 
the effective Lagrangian which describe the interaction between neutralino and quarks. If 
we consider only the extra Higgs exchange contribution, the effective Lagrangian reads 

Aff = A{{xx){ff) + Ai{xi'x){h'f) (3.22) 
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where 



^1 = ^^|^iViic^(Afi5S7 + ^16C7), (3.23) 
4ml 



A^ = -'^^^NnCj{Ni5S^-NiGC^), (3.24) 



■^i = -^^P^Nus^{Ni5S^ + iViec^), (3.25) 
4ml 

-42 = -^^?^A^llS7(iVl5S7-^16C7)- (3-26) 

4m^ 

Note that Ai^"^ comes from CP-even extra Higgs exchange and A2^^ comes from CP-odd 
extra Higgs exchange. Since the term with coefficient ^2^"^ §ive the velocity-dependent 
contribution to the elastic scattering, only the term with A^^'^ are important for our anal- 
ysis. 

The spin-independent cross section with a nucleus with Z protons and A — Z neutrons 
is given by 

4' = —{^pZ + Xn{A-Z)Y (3.27) 



where 



where Atv (-^ = P, n) is four-fermion coupling of the neutralino with point like nucleus and 
fj,^ = m-^mA/ {jnx ~^ "^-4) reduced neutralino-nucleus mass. Now we can obtain Xn 

from the parton level couplings as 

\n = {N\uu\N) At + {N\dd\N) Ai 

{N\ff\N) = "^ff, (3.29) 
m/ 

and we consider the contribution only from up-quark and down-quark because the contri- 
butions from other quarks are suppressed by their masses.^ The parameters and 
are obtained as [41] [42] 

/P = 0.023, = 0.033, = 0.018, /^^ = 0.042. (3.30) 

Since we can take arbitrary values for the Yukawa couplings of extra Higgs (y^ and 
y'^) independent of the masses of quarks, the effective coupling Ajv are generally different 
between proton and neutron. In other words, we can make large isospin violation for the 
spin-independent scattering as long as fj^ have different values each other [17] [18]. Now 
we parameterize the isospin violation by the ratio r = A„,/Ap and obtain 

y'd = k'^t^^'^i-y'u, k= hn [il - (3-31) 



^ The contribution from strange quark is possibly large if is much larger than and . But we 
adopt the small value of as suggested by the lattice study [43] . 
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Figure 1. Contour plot of y'^i/Uu'- We illustrate the contour for y'^/y'^ = ±[0.0, 0.4, 0.8, 1.2, 1.6] 
which is given by (3.31). 

In Figure 1, we show the contour plot of y'Jy'u for various values of r and 7 with the quark 
mass ratio rriu/rad ~ 0.55. Using (3.31), we can rewrite (3.28) as 

= ^iVii — - kf^)y'^c,{N^,s, + N^^c,). (3.32) 

To discuss the favored region for the direct detection, we often use the neutralino- 
proton cross section 

4 / my.m> 

^P=-( ^ 

vr \ -|- nip 

and the preferred regions of DAMA and CoGeNT signals are estimated to be 

1.5 X 10-^° cm2 < 4™^) < 3.0 X 10"^° cm^, (3.34) 

5.0 X 10-^1 cm2 < aJCoGcNT) < g_Q X ^^2 (3 35) 

for niy. ~ 8 GeV if there is no isospin violation (Ap = A„). Here we adopt the larger sodium 
quenching factor QNa = 0.5 it 0.1 so that the DAMA preferred region include the dark 
matter mass rriy. ~ 8 GeV [33] [34]. On the other hand, these parameter regions are already 
excluded by the experiments such as CDMS-II and XENONIOO. In particular, XENONIOO 
have reported a severe constraint 

CTp < 4.0 X 10"^^ cm^ (3.36) 



^'"^ ^ Xl (3.33) 



for ~ 8 GeV. To consider the case with an isospin violation, we introduce the effective 
cross section of neutralino and nucleus defined by 

a!k = V -^PiiXpZ + Xn{Ai - Z)f (3.37) 
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Figure 2. Scatter plot of direct detection cross section and experimental preferred/excluded 
regions for — 1 < r < 1. We set — 8 GeV, = 150 GeV, tan /3i = 5, = 1 TeV and analyze 
the parameter space 7 = [0, 27r), ^i2,/i2i = [100 GeV, 500 GeV) to draw the scatter plot. Yukawa 
couplings y'^ and y'^ are determined to generate the relic abundance of dark matter Q,^h^ ~ 0.1 
with the help of (3.19) and (3.31) In the scatter plot, we restrict |?/'J, \y'^\ < (47r)^/^. 

where i is summed over isotopes Ai with fractional number abundance pi. This definition 
reflects the practical constituents of the detector with assumption that only one of the 
elements dominates the contribution for the cross section. The effect of the isospin violation 
can be extracted by comparing with the cross section without isospin violation (r = 1): 



Therefore, the favored region of the neutralino-proton cross section in the case of isospin 
violation can be obtained by replacing Up with KrCTp. 

In Figure 2, we illustrate the scatter plot of direct detection cross section by varying the 
ratio of neutron to proton coupling r. Once we determine the ratio y'^/y'u from (3.31), the 
Yukawa couplings y'^ and y'^ can be derived by requiring to generate the relic abundance 
n^h? ~ 0.1 as computed in previous subsection. For evaluating dp, we adopt the mass 
of up-quark niu ~ 2MeV. We also show the preferred region for DAMA and CoGeNT 
and constraint from XENONIOO in figure 2. As can be seen from the figure, this model 
can generate the cross section consistent with the overlap region of DAMA and CoGeNT 
experiments. Moreover, the constraint from XENONIOO is weakened for r ~ —0.7 and we 
can find the parameter region consistent with these three experiments. We have to mention 
that CDMS-II experiment still exclude the region favored by DAMA and CoGeNT even if 
we introduce isospin violation. This contradiction should be investigated carefully in the 
future experiments. Recently, CRESST experiment has reported that they found signal 





(3.38) 
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consistent with the light dark matter scenario [3]. We can also see that the 3a preferred 
region for CRESST result is compatible with that of DAMA and CoGeNT if we take an 
isospin violation r ~ —0.7. 

4 Collider signature 

4.1 Constraints from precision experiments 

Although we have managed to explain the results from dark matter experiments by su- 
persymmetric four-Higgs doublet model, there can be some constraints from the precision 
experiments for the parameters in this model. In this subsection, we give some comments 
on the possible constraints from the current experiments. 

Since we now consider a very light neutralino with mass ~ 8 GeV, we should care 
about the constraint from the invisible decay rate of Z-boson. The decay rate of Z-boson 
to a pair of the lightest neutralino can be written as 



which is similar to the MSSM case [44] [45]. Therefore, arbitrary values of A'^is and Niq 
given in (3.11) and (3.12) with |/ii2|, |Ai2i| ^ 100 GeV are consistent with the constraint 
r(Z° ^ XiXi) < 3MeV as discussed in [21]. 

When we introduce the extra Yukawa couplings, the severe constraints from flavor 
physics generally appear as pointed out in the previous section. Although the tree-level 
contribution for the flavor changing process is absent by setting the Yukawa couplings of 
extra Higgs fields to the form given in (3.15), we should also care about the loop contribu- 
tions mediated by extra Higgs fields. In this paper, we simply assume that the dangerous 
processes such as 6 — t- s 7 are suppressed by the suitable choice of the off-diagonal elements 
of extra Yukawa couplings and supersymmetry breaking parameters. 

4.2 Wjj anomaly 

As discussed in the previous section, this model predicts light neutral extra Higgs fields 
(j)^^ and cpQ^ to account for the relic abundance and the preferred region for the direct 
detection of dark matter. If such light fields exist with large Yukawa couplings, there is 
a possibility that the current experiments already have an ability to discover (or exclude) 
these particles. On the other hand, the CDF collaboration has reported that there is an 
excess for the dijet mass distribution of the W + j j event around 150 GeV. This excess 
may be explained in our model if the final state two-jets are mediated by the extra Higgs 
field. Actually similar situation has been investigated in [37] which introduce an additional 
quasi-inert Higgs doublet to the standard model. Although this anomaly may not be 
true because no anomaly is seen by DO experiment, it is still important to investigate the 
possibility to test the model by the existing experiments. 

When we consider the dark matter relic abundance and the rate of direct detection, 
the coupling constants are evaluated around the mass scale of the neutralino m^. Now we 




(4.1) 



- 13 - 



2.0 



0.5 



[GeVj 



Figure 3. Cross section ctq defined by (4.5) as a function of the mass of extra neutral Higgs m^. 

have to incorporate the effect of the renormahzation group flow to calculate the production 
cross section of the extra Higgs. Below the mass scale of the extra Higgs, only the SU{3)c 
gauge coupling gives significant contribution for the running of the Yukawa couplings y'^ 
and y'^. So the renormahzation group equations we should consider are given as 

at one-loop level. ^ These equations can be solved to obtain 

• y'u/di^x) (4-4) 

where Aa = as{m).) — as{m^) with as = g1/{AT:). For m;^ = 8GeV and = 150 GeV, 
we obtain ^^/^^(w-t/)) — 0.7 • Uu/dimx) ^om as{m^) ~ 0.2 and as(m0) ~ 0.1. 

If we accept the result from the CDF collaboration as the true signal of new particle, 
their result suggests that the size of new physics contribution for the cross section of 
W + j j process is comparable with that from W'^ W~ production process. Since a{pp — )• 

W^) ~ 8.6 pb for ^/s = 1.96 TeV, it roughly require the cross section for the process 
pp — )• (pEi/oi around a few pb. 

To estimate the cross section, we assume that extra charged Higgs are much heavier 
than and consider only the f-channel quark exchange contribution of (j)Ei and 4>oi final 

^ More precisely, we consider the running of the effective couplings A{{xx){f.f) + -^2 (x7^x)(/7^/) ^-nd 
take a matching at the scale m^. 



y'u/di.^<f>) 



l + O 



4tt 
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states for simplicity. Note that the contribution from the resonant production of charged 
extra Higgs (j)^ is not so large because the branching ratio Br(0^ — )• (fP) is very small 
compared with Br((^^ — t- j j) through large Yukawa coupling. We write the cross section 
of this process as 

'y{pp ^ W+ ^ [y'^a, + y^^) (4.5) 

and (To for various mass of extra Higgs is plotted in Figure 3. Here we used LHAPDF 
[46] to compute the cross section. From this figure, we can see that required size of cross 
section can be obtained by 0(1) Yukawa couplings. Of course, we have to carry out a 
Monte Carlo simulation to check whether this model can generate a signal consistent with 
the experiment, but it is beyond the scope of this paper. Since this model also gives large 
cross section for the processes pp ^ (pEi/oi and pp ^ j 't'Ei/oi^ precise measurements 
for these processes are important to test the model. 

5 Conclusion 

In this paper, we have considered supersymmetric four-Higgs doublet model to explain the 
signals of dark matter direct detection from DAMA and CoGeNT experiments. We have 
evaluated the relic abundance and dark matter-nucleus cross section for the neutralino LSP 
scenario including the contributions from the extra Higgs superfields. 

We have shown that the preferred regions given by DAMA and CoGeNT can be simul- 
taneously explained by the supersymmetric four-Higgs doublet model with the parameters 
which are consistent with the observed value of dark matter relic abundance. To explain 
the results of DAMA and CoGeNT, an isospin violation for the neutralino- nucleus scatter- 
ing is essential and we can achieve this situation by adjusting the ratio of extra Yukawa 
couplings y'Jy'u- 

The extra Higgs fields we introduced as an explanation of the light dark matter are 
candidate for the solution of Wjj anomaly reported by CDF. We have calculated the cross 
section for the process pp ^ 4>^ei/oi ^^'^ ^^^^ ^^^^ cross section is compatible 
with the experimental data. 
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